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Explosive particle production due to parametric resonance is a crucial feature of reheating in infla- 
tionary cosmology. Coherent oscillations of the inflaton field act as a periodically varying mass in 
the evolution equation for matter fields which couple to the inflaton. This in turn results in the para- 
metric resonance instability. Thermal and quantum noise will lead to a nonperiodic perturbation in 
the mass. We study the resulting equation for the evolution of matter fields and demonstrate that 
noise (at least if it is temporally uncorrelated) will increase the rate of particle production. We also 
estimate the limits on the magnitude of the noise for which the resonant behavior is qualitatively 
unchanged. 
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I. INTRODUCTION 



Over the past few years, it has been realized that explosive particle production at the end of the period of inflation 
is a crucial aspect of inflationary cosmology. At the end of the period of exponential expansion, the energy density of 
matter and radiation is exponentially small, and without a very fast transfer of energy from the inflaton to ordinary 
matter, it is not possible to obtain a high post-inflationary temperature. 

As was first pointed out in B and discussed in more detail in |^]-^| and many other recent papers @ , an oscillating 
scalar field induces a parametric resonance instability in the mode equations of any bosonic matter fields which couple 
to it. In particular, this applies to the inflaton, the scalar field responsible for inflation. At the end of the period of 
exponential expansion of the Universe, the inflaton is predicted to be performing homogeneous oscillations about its 
vacuum state. This will lead to instabilities in the mode equations of any bosonic matter field which couples to the 
inflaton, and this instability corresponds to explosive particle production. 

Instabilities, due to resonance effects, are in general quite sensitive to the presence or absence of noise. In a 
real physical system we expect some non-periodic noise in the evolution of the inflaton. Such noise could be due to 
quantum fluctuations (the same quantum zero point oscillations which are believed to be the source of classical density 
perturbations in the inflationary Universe scenario) or thermal fluctuations. [] In both cases, the amplitude of this 
noise is expected to be very small. Nevertheless, it is important to analyse the effects of this noise on parametric 
resonance. 

In this paper, we study the effects of noise in the inflaton field on the evolution equation of matter fields which 
couple to the inflaton. As a first step, we take the noise to be homogeneous in space. In a subsequent paper we 
plan to study the (more realistic) case of inhomogeneous noise. Our main result is that noise which is not correlated 
temporally will on the average lead to an increase in the rate of particle production. We also derive limits on the 
amplitude of the noise for which the parametric resonance behavior persists. This result eliminates a further doubt on 
the effectiveness of resonance in realistic inflationary models. Most of our results apply for both narrow and broad- 
band resonance. For simplicity we neglect the expansion of the Universe; however, we do not believe that including 
effects of expansion would change our main conclusions. 

1 Note that in the usual analysis, the inflaton is treated as a classical background field. 
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In this paper we shall consider a simple model of reheating in the presence of a spatially homogeneous noise term. 
The inflaton field is taken to be a real scalar field and is denoted by tp. In the period immediately after inflation, 
<p(t) is assumed to be oscillating coherently about (one of) its ground state(s). We shall, however, include a small 
aperiodic perturbation, i.e. 

<p(t) = Acos{ut) + q(t) , (1) 

where lo is the natural frequency of oscillation of tp and q(t) denotes the noise. 

We shall (following Jl|-f| ) take p to be coupled to a second scalar field \ (which represents matter) via an interaction 
Lagrangian which is quadratic in \, for example of the form 

C-%nt oc ^px 2 ■ (2) 

To simplify the discussion, we neglect nonlinearities in the equation of motion for x- Such nonlinearities may very 
well be important and lead to an early termination of parametric resonance. This issue has recently been discussed 
extensively, see e.g. but is not the topic of our paper. Instead, we are interested whether the presence of noise 

such as included in ([!]) will effect the onset of resonance. 

For simplicity, we will neglect the expansion of the Universe. In it has been shown that the expansion can 

be included without difficulty, and that it does not prevent the onset of the parametric resonance instability. Since 
the equation of motion for x{ x i t) is linear and translation invariant, the Fourier modes evolve independently. We can 
immediately write down the evolution equations for the Fourier modes of denoted by Xk, in the case of "coherent" 
noise. They take the form 

Xk + H + (p(wi) + q(t))] Xk = 0, (3) 

where p(x) is a periodic function whose frequency is 2tt, q(t) represents the noise which we consider as a perturbation 
of the driving function p(t) , and 

u 2 k =m 2 x + k 2 . (4) 

Note that our analysis applies both to the case of narrow-band fli (p/uj 2 -C 1) and broad-band @ (p/uj 2 ^> 1) 
resonance. 

By assumption, q(t) is a small perturbation of p(t). Hence, it is suppressed compared to p(t), and in order to 
represent this, we write 



q(t) = guj 2 n(t) , 



(5) 



where n(t) is a dimensionless "noise function" of amplitude 1, and where we have introduced the dimensionless small 
coupling g by extracting the dimensions of q by means of the factor u> 2 . 

We wish to consider how the presence of noise in the inflaton effects the excitation of the matter field components 
Xk from their vacuum state. Another way in which thermal noise can effect the calculation has recently been studied 
in where it was shown that parametric resonance is also effective if we consider the field \ to be initially in a 
thermal state. 

We can study this problem in three ways (see e.g. for a recent review). First, we study the effects of the 
noise on the solutions of the classical equation of motion (||) using the method of successive approximations and the 
Furstenberg theorem on products of random matrices. Next, we use the Bogoliubov mode mixing technique @ to 
analyze the change in the solutions due to the noise. This method is closely related to a consistent semiclassical 
analysis which treats the excitation of the x field as a problem of quantum field theory in a classical background 
inflaton field (see e.g. Appendix B of ||). Finally, we study the effects of the noise using the Born approximation 
and compare with numerical results. 



The principal results of this section are that the exponential growth rate of solutions is a continuous quantity 
with respect to q(t), so that a small addition of noise will not change it overly. Furthermore, under reasonable 
assumptions of decorrelation of the noise, the growth rate is shown to be always strictly increased by the noise. For 
these considerations it proves convenient to rephrase our basic second order differential equation (|) as a first order 
2x2 matrix differential equation 



II. EXACT RESULTS 



$ 9 = M(q(t),t)$ q 



(6) 



with initial conditions $ 9 (0, 0) = /. Here, M(q(t),t) is the matrix 




(7) 



and ^(t, a) is the fundamental solution (or transfer) matrix from time a to time r; 




(8) 
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consisting of two independent solutions 4>i(t;q) and </>2{t;q) of the second order equation (|3|). 

For vanishing noise, i.e. q(t) — 0, the content of Floquet theory is that the solution of @ can be written in the 
form 

*o(*,0) - Po(t)e ct , (9) 

where Po(t) is a periodic matrix function with period T = and C is a constant matrix whose spectrum in a 

resonance region is spec(C) — {±/i(0)}. 

We would like the noise q(t) to correspond to random quantum or thermal fluctuations superimposed on the 
periodic classical oscillation of the inflaton field if(t). Our picture is that of fluctuations driven by Brownian motion, 
however for the purposes of deriving the exact results to be presented here, it is sufficient to make certain statistical 
assumptions about the noise q(t). These are phrased in terms of a sample space 51 from which the realizations k of 
the noise q(t) — q(t; k) are drawn. On the sample space f2 there is a probability measure cLP(k), and expectation 
values of functions f(n) with respect to this measure are denoted by 

E(f) = J dP{n)f{K) . 

For our purposes we may take = C(R), the space of bounded continuous functions on R, and cLP(k) on Q a 
translation invariant measure. We assume that the noise is ergodic, which is to say that 
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E(f(q)) = lim - / dr /(?(•+ t;/s )) (10) 







for almost all dP(n) realizations ko of the noise. 

The first relevant mathematical result (Proposition 1, see Appendix) is that the growth rate (the generalized Floquet 
exponent, or Lyapunov exponent) of the solutions of (^) is well defined by the limit 

M(g) = Jim ^ log ||nf =1 $(jT, (j - l)T)\\ , (11) 

where || • || denotes some matrix norm (the dependence on the specific norm drops out in the large N limit). The limit 



in (11) exists for almost every sample (with respect to dP(n)), and it is almost everywhere dP(n) constant, which is 
to say that it depends only upon the the statistics of samples q(t] k) and not on the individual realizations. 

The first qualitative result (Proposition 2, see Appendix) is that the growth rate n(q) is continuous in q in an 
appropriate topology on Q (the topology on C(R) of uniform convergence on compact sets). In particular, this 



implies that if we write the noise as in (|5|) in terms of a dimensionless strength g and view the statistics of the 
noise as a measure on the function space of n(t), then /i(q) converges to /i(0) as g converges to zero. We remark 
that Propositions 1 and 2 do not depend on the fact that the system is one-dimensional, and will thus also hold for 
inhomogeneous noise. 



The limit (11) also describes /i(0), which is implicitly a function of the background periodic potential p(t); it shows 
in particular that /i(0) > 0. A second qualitative estimate (Theorem 3, see Appendix and jl8|) is that, if the support 
of the probability measure dP(n) on the space C(R) includes the sample q(t) = 0, then ji(q) > whenever /Lt(0) > 0. 

The results cited above do not give insight into the quantitative value of the growth rate in the presence of noise. 
A first step in this direction is given by the following theorem. Suppose that the statistics of the noise are such that 

(i) The noise q(t) is uncorrelated in time on scales larger than T, that is, {q(t; n) : jT < t < (j + 1)T} is independent 

of {q(t; k) : IT < t < (I + 1)T} for integers I ^ j, and is identically distributed. 

(ii) Restricting the noise q(t; k) to the time interval < t < T, the samples {q(t; n) : < t < T} within the support 

of the probability measure fill a neighborhood, in C(0, T), of the origin. 

Hypothesis (i) implies that the noise is ergodic, and therefore the generalized Floquet exponent is well defined. We 
will show (Theorem 4, Appendix) that in fact fi(q) is strictly larger than /x(0); 

//(<?)> M(0), (12) 

which demonstrates that the presence of noise leads to a strict increase in the rate of particle production. This is a 
quantitative lower bound on the generalized Floquet exponent with noise. Our result is based on an application of 
Furstenberg's theorem, concerning the Lyapunov exponent of products of independent identically distributed random 
matrices {tyj : j = 1, TV}. It states that (modulo certain assumptions which are shown to hold in the Appendix) 

j^lin?=i<M = A>o, (is) 

where A = X(q) depends again only on the statistics of ^j, and not on the individual samples. A further result is that 
for any nonzero vectors v\ and V2 and for almost all 

lim J-logKlFjLi*^) = A. (14) 

N— too iv 1 



In order to apply Furstenberg's theorem to obtain (|12|), we start by factoring out from the transfer matrix $ g (t,0) 
the contribution due to the evolution without noise; 



*,(i,0) = * (*,0)*,(t,0) = P (i)e c %M) • (15) 



The reduced transfer matrix ^ q (t,0) satisfies the following equation 



-q{t;n) 



*, = S(t; K )V q = $^(i) ( _ ° , ° ) (16) 



which can be written as a matrix integral equation 



%(t,0) =1+1 dr 5(r;/c)*,(r,0). (17) 



Solutions ^ q {t, 0) to (17) are constructed in the next section, using the method of successive approximation, giving 
rise to the transfer matrices ^ q {jT, (J — 1)T) which are the fundamental solution matrices for the period intervals 
of the background potential. By properties (ii) of decorrelation of the noise, the quantities ^f q (jT, (j — 1)T) are 
independent and identically distributed for different integers j, and we can apply the Furstenberg theorem to the 
following decomposition of *f? q (NT, 0); 

V q (NT,0) = nf =1 * ? (jT, (j - 1)T) . 

We may choose for instance the vector Vi in (Q) to be an eigenvector of $o(T, 0)' = (Po(T)e CT y , the transpose of 
the transfer matrix of the system without noise, with eigenvalue e >1 ^ T . Then dl4| ) becomes 

^ log 1(^,^(^,0)^)1 = J Llog\(v 1 ,Po(NT)e CNT * q (NT,0)v 2 )\ 

= ^log(e^ JVT |( Wl ,vI/ 9 (7VT,0> 2 )|) (18) 
= M0) + ^log|( Vl ,nf =1 vI/ g (jT,(j-l)T)« 2 )|. 

Taking the limit TV — > oo and applying ( |l3| ) and (|lj) we obtain 

fl(q) = fx(0) + A > fl(0) (19) 

which proves the main result (12). Note in particular that ( |l2] ) implies that even for modes which without noise are 
in a stability band (/i(0) = 0) there will be exponential particle production in the presence of noise. The fact that 
there is particle production due to the presence of noise is a nontrivial result, and the fact that the rate of particle 
production is exponential may be understood physically as a sign of stimulated emission. Note that, strictly speaking, 
the ergodic hypothesis is only satisfied in the t — > oo limit. Finite time intervals, such as the characteristic time of 
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reheating /U -1 , may be insufficient for the decorrelation of the noise. In this case, we may observe transient effects 
(see Section VII). 

So far, our results have not given any quantitative upper bound for the growth rate of modes of x m the presence 
of noise. Additionally, they do not show how the evolution of an individual solution is modified, as compared to what 
occurs in the system without noise. In the following sections we will give estimates on the magnitude of noise for 
which we can demonstrate that the parametric resonance behavior remains qualitatively unchanged. These estimates 
are then used in the proof in the appendix of the main result ([[I). 

III. AN ESTIMATE USING SUCCESSIVE APPROXIMATIONS 



In this section we sketch the construction of the fundamental solution matrix *5f q (t,0) for system (16), showing in 
particular that {^ q (jT, (j — 1)T)}°^ 1 represents a sequence of independent identically distributed random matrices. 
The starting point for this estimate is the integral equation ( )l7| ) for the reduced transfer matrix $(i, 0) (dropping the 
subscript q for notational ease). We solve this equation by constructing a sequence of approximate solutions 

*„+i(t,0) = /+ / dr 5(t;?(t))*„(t,0) (20) 
Jo 

with = I- The differences between successive terms of the sequence satisfy 

* n+ i(i,0)-* n (i,0) = / dTS(T;g(T))(* n (T,0)-* n _i(T,0)). (21) 



t 



We will apply this equation for evolution over the time period T. 

By induction is can easily be shown that the successive differences satisfy the estimate 

||*n+i(r,0) - *„(T,0)|| < M n+1 — — , (22) 

[n+iy. 

where the constant M is 

M = sup ||S(s;<?)||. (23) 

0<s<T 

From the definition of the matrix S (see jl^ ) and (||)) it follows that (if the norm || ■ || is taken to be the supremum 
norm or the norm mentioned after (|ll|)) 

M < e 2 ^ 0)T \\q{--K)\\. (24) 
The following telescoping sum then describes the fundamental solution matrix; 



*(t,0) = 5^(* n (t,0)-* n _i(t,0))+/, (25) 



and by using (^) we estimate that 

(Aft)™ 



n*(*,o)-/|i < 



n 
n=l 



= e M * - 1 < (Mi)e Mt . (26) 

The fundamental solution of (||) is $ q (t, 0) = $o(^ 0)^(4, 0), therefore the deviation between the solution with noise 
and the solution without it can be measured by the estimate 

||3 9 (i,0)-*o(t,0)|| = p (t,0)(* 9 (t,0)-i)[| (27) 
< ||$o(t,0)||||iJ s (t,0)|| , 

where R q (t, 0) = (&q(t, 0) — I). When we set q(t; k) — gu> 2 n{t; k), and g is taken to zero, then the constant M in ( Kj ) 
also converges to zero, therefore for fixed time t the error 0)|| < (Mt)exp(Mt) is small and thus solutions of 

the two equations which have the same initial data are shown to be close. 

We will now make use of the above estimate in order to bound from above the difference n(q) — /x(0). We consider 
evolution over a time interval NT and break up the transfer matrix into matrices corresponding to single periods T. 
From (P) and (jll]) we obtain 



tiQ)-m = ^^(logH ||nf =ie TC|| ]) • ( 28 ) 

Making use of the periodicity of Po, we can write the norm of the product in the numerator (denoted by NUM) as 

NUM = \\uf =1 (e TC + A q (jT,(j- 1)T))\\ , 

where A q represents the difference in between the transfer matrices with and without noise. We need to bound its 
contribution in magnitude from above. Due to the fact that the noise is uncorrelated over time intervals larger than 
T, this needs to be done only for one period. Introducing the symbol 



we can bound the difference in the generalized Floquet exponents (|2q) by 



M?) ~ M(0) < Jim -±= log(l + 9) w = I log(l + 6) , (30) 



which for low amplitude noise can be approximated by 0/T. 



Making use of the fact that /x(0)/w <C 1, it follows that for noise of the form given by (||), Q is of the order gui 2 T 2 ^ 
Hence, from ([30]) we conclude that if 
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M(0) 



(31) 



the effects of the noise do not significantly effect parametric resonance. Furthermore, the estimate ( p2| ) shows that 
solutions of (6) cannot deviate by more than MT exp(MT) from §o(t, 0), over time intervals of length T. 

IV. A SECOND ESTIMATE USING SUCCESSIVE APPROXIMATIONS 

The growth rate fj,(q) of the previous section is the result of a limiting process over many periods T of the background 
periodic excitation pit). It may be, however, that the noise acts over small characteristic subintervals of length At, and 
the reheating epoch only encompasses a finite and relatively small number of periods of the background. By studying 
directly the second order differential equation (3) with an estimate based on the method of successive approximations, 
we will exhibit in such a situation a bound in terms of the coupling constant g on the deviation of solutions perturbed 
by the presence of noise from the unperturbed solution. In case this deviation is small when compared with the size 
of the unperturbed solution, we can conclude that the growth due to parametric resonance is not destroyed by the 
presence of the noise term. We obtain an estimate on g for which this holds; it turns out to be of the same character 
as estimate (|3l|). 

The method is to rewrite equation (^) as an integral equation and to determine the general solution in terms of a 
series of successive approximations. In the following we use a well-known result for differential equations. Let x{t) be 
the solution of the differential equation (we drop the index k on x to simplify the notation) 



2 One factor of T comes from the time integration in the transfer matrix, the second arises since a factor of T must be inserted 
in the matrix element in (jig) if both components of ^ are to have the same dimensions. 



X+H + (P(t) + q(t))] X = 



(32) 



subject to given initial conditions, and let Xh be the solution of the unperturbed equation 



X + [wfc+p(t)]x = 0, 



(33) 











satisfying the same initial conditions. Then, Eq. ( |32| ) can be rewritten as an inhomogeneous Volterra equation liS of 
the second kind 



x(t)=Xh(t) + 



w 



-q(t )x(t 



dt , 



(34) 



where <pi and (f> 2 are two independent solutions of equation ( p3| ) and W is their Wronskian, which is independent of 
time (by Abell's formula, see e.g. ]15|]). We then apply the method of successive approximations starting with the 
unperturbed solution Xh(t). We define recursively 



Xo(t) = Xh(t) , 

Xn(t)=Xh(t) + 



w 



-q{t )Xn-l(t ) 



dt n > 1 



(35) 



From the above, we get an appropriate estimate for the successive approximations 



Xn+l{t) ~ Xn[t) 



w 



-?(*) Xn{t ) -Xn-l{t ) 



dt . 



(36) 



The function p(t) is periodic, so that the two independent solutions of equation (|33J) can be written in Floquet's 
form @ 



Mt) = e"<°>*pi(t) , fc(t) = e-^ t P2 {t) , 



(37) 



where pi and P2 are both periodic functions of time with period T, and /i(0) is the Floquet exponent discussed in the 
previous section. Let us focus on the behavior of the solution over an arbitrary time interval I = [ti,tt], of length 
AT. Define 



Mi = max \pi (t) I , M 2 = max \p 2 (t) 

0<t<T 0<t<T 



According to our assumptions, the perturbation q(t) is bounded by guj 2 , 



max \q(t)\ < guj < +oo . 
t gr 



(38) 



(39) 



Defining the norm 



\Xn - Xn-i|| = max|x„(t) - Xn-i(t)\ , 



(40) 



and using equations (|36|)-(p9|) it follows that for t £ I then 
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- x.MI < ( 2 .'" 2 ^<°>-)" +1 <i^| ta || • <«) 

Define 

a = 2 5W 2 ^e^. (42) 
Let us now estimate the distance d — ||x — Xftll- We consider the telescopic series 

n 

Xn+i(t) - Xh{t) = -Xj) > (43) 

4=0 

and use equations (f4l| ) and (ji|) and the triangle inequality to obtain 

™ " (t — t V+ 1 

\ Xn+1 (t)- Xh (t)\ < *£\ Xj+1 ( t )- Xj (t)\ < J2 a ' + c + lV WXhW, (44) 

3=0 j=0 V + 



where a is given by (|42j). Taking the limit n — > oo and using the definition of the norm in ( |40| ) we obtain 

llx - Xftll < (e QAT - l)||Xfc|| < aATe aAT \\ Xh \\ . (45) 

The last step is to ensure that the unperturbed solution has undergone appreciable growth by the end of the 
considered time interval. Typical initial data for the problem is of the form 

Xh(t) = CiMt) + (46) 

with coefficients c\ and C2 which are of the order 1, therefore = O(exp(/i(0)AT)). 

Note that the period of reheating in the absence of the noise is 7V/i(0) -1 , with N being a number of the order 
2 x \og{r)/H) (where H is the expansion rate of the Universe and ij is the energy scale of inflation). A typical value 
of T]/H is of the order 10 6 , and this implies that the number iV is not necessarily large. Hence, for practical purposes 
it is sufficient to control the effects of the noise over a time interval ^(O) -1 . In addition, in order to prove that 
the effects of the noise are small compared to the unperturbed solution, we must consider time intervals over which 
the unperturbed solution has significant exponential growth, i.e. over which the growing mode solution of equation 
(|33|), 4>i(t), dominates the decaying mode (j>2(t). This provides another reason to consider time intervals such that 
H(0)AT > 1. To be definite we will take ^(0)AT = 2. 

The maximal distance d between X and X h is controlled by ( |45| ) over any interval /. If the amplitude of the noise 
is sufficiently small, then this distance will not be appreciable when compared with the size of X h- In order for this 
to be true, we require that 
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aATe aAT < ~ . (47) 

The condition of smallness which emerges from ( p7|) is that 

a AT = 2e 2 gLu 2 (M 1 M 2 /\W\)AT « 1 . (48) 

To obtain an estimate of the order of magnitude of this constraint, we note that \W\ <~ MiM^uj. Inserting this into 
( p8| ) we obtain 

9«^, (49) 
w 

which is the same requirement as was obtained in the previous section (|3l|). Thus, in this case the parametric resonance 
growth of Xh(t) is preserved in x(t). 

Certainly there are better estimates that can be obtained if more information about the noise is provided. On 
the other hand, we have shown that for any kind of small homogeneous noise, over time intervals of length AT = 
O(/i(0) _1 ), parametric resonance is not destroyed by the presence of low amplitude noise in the inflaton field. 

V. AN ANALYSIS VIA THE BOGOLIUBOV METHOD 

Another approach which can be used to demonstrate that random noise with a small amplitude does not eliminate 
the parametric resonance instability is the Bogoliubov method @ . This approach has the advantage of giving some 
information on how the noise effects the dynamics. However, it is less rigorous than the method of successive 
approximations and holds only in the case of narrow resonance regime. 

The evolution equation for Xk{t) in the presence of the oscillating inflaton field (see Eq. (|3|)) is (dropping the 
subscript k and introducing an explicit expansion parameter e <C 1) 

X(t) + [u 2 k + e(p(ut) + <?(*))] X (t) = , (50) 

where p(u>t) is a function with period 2tt/lu and q(t) represents the noise. 

For simplicity we consider a mode k which (in the absence of noise) is in the first resonance band, i.e. for which 
= uj — A, where A is a small quantity. In this case we may write: 

p(ut) = ujI cos(2w fc + A)t , (51) 
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with — euik < 2A < eui k - Following the Bogoliubov method we assume that, to first order in e, the solution of 
( |50| ) is of the form 

A A 

X (t) = a(t) cos(w fe + —)t + b(t) sm(uj k + —)t , (52) 

where a(t) and b(t) vary slowly in comparison with the trigonometric factors. Notice that the exact solution includes 
terms with frequencies higher than uj k + A/2. These terms are neglected in the first order solution since they are of 
higher order in e. 

Inserting ( |52j ) and (|5l|) in ( ^o|) and keeping only first order terms we obtain 

; A e o \ . A, 



+ ( b — Loa — —tub — -0J^.b + eqb ) sin(a;fc + —)t = . (53) 



—tub - -u k b + eqb I sm(w fe + — 

It is clear that the coefficients of both the cosine and sine in the above equation must be zero, which yields two 
differential equations for the functions a(t) and b(t). We seek solutions of the form 

a(t) = A(t)e"* , b{t) = B{t)e^ , (54) 

where fi is a constant of order e. We then get 

.... As 
A + 2fiA + ujB + [iuB - ~^uA + -w\A + eqA = , 

.... A e 

B + 2[iB -ujA- ufiA - —ujB - -u 2 k B + eqB = . (55) 

Let us initially study the case when the noise is neglected. Since q — 0, by following the Bogoliubov method we 
can assume that A, B, A and B are of second order in e. The resulting system is (recalling that u> = 2ojk + A) 

2{j,A + BA + ^uj k B = 0, 2fiB - AA + ^w k A = , (56) 

which yields 

. ,2 / jr 2 A2 

2 _uj k I e A 



11 4 [l ~ 14 ) ' A ~ 5 = C °" St (57) 

This result means that as long as /x is real (— euj k < 2A < euj k ) and for g = the Xfc field grows exponentially in time. 

From the above analysis a simple but very important result follows. In the narrow resonance regime, the resonance 
is not destroyed by noise which is small compared to the amplitude of the oscillatory inflaton field p(u)t). Specifically, 
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if g (see (|])) is small, say g = 0(e), then the last term in both of the equations ( |55| ) can be neglected and the result 
of Eqs. d57|) follows. 

On the other hand, when q is not negligible, then the derivatives A, B 7 A and B must be of first order in e and 
not 0(e 2 ) as we assumed above. This can be seen as follows. Solving equations ( p5| ) as before under the assumption 
that the derivatives are negligible leads to a solution for A and B which depends on time (via the noise q), and 
whose derivatives are thus proportional to eT, where T is the rate of change of the noise. This demonstrates that the 
derivatives are not negligible. 



In order to circumvent this difficulty we differentiate once equations (55) and neglect all terms of second and higher 
order in e to get the following 

A + loB + eqA = , 

B - ujA + eqB = . (58) 
By integrating the above system we get (to first order in e) 

A + uoB + eqA = Ci , 

B-LuA + eqB = C 2l (59) 

where C\ and Ci are both constants of order of e. 

The next step is to check the compatibility between (|59| ) and (|5q). Inserting ( |59"|) into ( [55| ) and neglecting terms 
of order e 2 , we obtain a matrix equation for the vector (A,B) which has a solution only if the determinant of the 
coefficient matrix vanishes. This condition leads to the same value of the Floquet exponent /j, as obtained in the 
absence of noise (cf. Eq. (|57|)). Since q(t) is a general time dependent function, a solution exists only if C\ = Ci = 0, 
and only if the terms ujB and w A in ( |59| ) are negligible compared to the other terms. Any other choice for C\ and C*2 
implies A = const., B = const., which is only the case for g = 0(e) as we showed above. 



Moreover, using ( pq ) and (59) it is seen that conditions \toB\ <C \sqA\ and \ujA\ <C \eqB\ are equivalent to w C 
|(?|/|(7| ~ r, where T is the characteristic rate of the noise. In this situation, Equations ( |59"| ) reduce to A + eqA = , 
B + eqB = , whose solutions to first order in e can be written as A , B cx e^ 1 ^, where j(t) = e J ^(r)rfr with 
* = q = guj 2 n(t) (see Eq. (|)). 

In this case, parametric resonance can be destroyed by possible exponential decay of A and B. For assume that t 
is sufficiently big and that the noise n(t) is a random walk, so that we can define the average of n in the time interval 
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from to t as n = \ f o n(r)dr - l/^r) 1 / 2 , which yields j(t) ~ geu 2 (t 3 /T) 1 ' 2 . This leads to A, B oc e -^ 2 «* 3/2 r 1/2 ; 
which implies that (see Eqs. a(i) and b(t) vary exponentially with /it — j(t). Since 7(t) grows faster than t, for 

sufficiently large times and for n > it follows /it — j(t) < 0. To be precise, let us assume that j(t) is positive and 
define the time t e such that pd, e — 7(i e ) = 0. This furnishes t e ~ r/g 2 cu 2 . Therefore, for t > t e , a(t) and b(t) decay 
exponentially with time, destroying the resonance. However, for n < the resonance is enhanced by the noise. 

Notice that eui is of the same order as fj, so that the noise is effective for times t > t e , where t e is such that 
g 2 u>t e = T/ui ^> 1. For g sufficiently small, t e can be bigger than the whole reheating time, implying no effect on 
the exponential growth of the resonant field x- This result is consistent with the results of the previous sections and 
also with the numerical studies reported in Sect. VII (see also section VI). Based on the solution without noise, the 
reheating time can be taken to be of the order fi . Thus, a sufficient condition for noise not to prevent the onset of 



resonance is 

5 2 «^-. (60) 

UJ UJ 

Since T > /i, this condition is less restrictive than the result (fl9]) obtained in the previous section. Note, in particular, 

that the larger T is, the less sensitive the resonance is to the effects of the noise. 

Concerning the opposite limit u> 3> l^l/l^l to which the above analysis is not applicable, there is however an alternate 

procedure to study the problem. In such a case, the periodic function oscillates many times during a characteristic 

time step of the noise T . The noise can then be thought as a slowly varying function when compared to p{ujt). 

Therefore, q(t) is "almost" constant within a time interval of the order of T -1 and can be considered as a perturbation 

on the frequency uj 2 . The analysis via Bogoliubov method can be repeated for each of the N = tT time intervals, but 

with a different frequency squared oj 2 + eqj for each interval. Here qj can be chosen as the extremum value of the 

1/22 / \ 2 s 



noise within the j-th time interval. The Floquet exponent is then /1 — — — I A + s— . We see that the 

4 V 4 V UkJ J 

resonance is preserved provided that the inequality —ELUk < 2(A + sqj/uJk) < su>k holds for all time intervals. The 
average Floquet exponent is obtained by replacing qj in the above expression by its mean value over all reheating 
time q ~ gu\j '(tT) 1 / 2 . Notice that (for fixed uik) \<j\, and then the effect of the noise, decreases with V and t. This is 
also verified numerically (see Sect. VII) and it holds approximately even in the case ujk — w/2 <C |?|/M- 
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VI. AN ESTIMATE BY MEANS OF THE BORN APPROXIMATION 



The two previous methods give (estimates for) lower bounds on the amplitude of the noise in the infiaton field 
for which we can prove that parametric resonance persists. However, it is to be expected that resonance persists for 
substantially larger amplitudes. In this section, we adopt a perturbative technique to estimate the strength of the 
noise required to change the resonant behavior of the modes. 

The starting point is the mode equation (^). We will solve this equation in the first order Born approximation. We 
write the solution \ (dropping the mode index k) as 

X = Xh+Xq, (61) 

where Xh is the solution without noise satisfying the given initial conditions, and Xq is the contribution of the noise 
to x computed to first order in the Born approximation, i.e. satisfying the equation 

Xq + {u% + ep{t)) Xq = eq{t) Xh (t) (62) 

and with vanishing initial data. If we introduce a noise coefficient g as in (^|), then our approximation corresponds to 
first order perturbation theory in g. 

For the time dependence of ([!]), i.e. p(t) ~ cos(ut), the "homogeneous" solution Xh can be (to first order in e) 
written as 

Xh(t) = cie"*cos(|i + ipx) + c 2 e-^cos(^t + ip 2 ) , (63) 

where c\ and c 2 are the amplitudes of the two fundamental solutions of the homogeneous equation (denoted Xi W 
and X2(t)), and ipx and (p 2 are phases. 

By means of the Greens function method, the solution of (|62|) takes the form 

X q (t) - Xi(t) f dt'W{t l r 1 X2(t'W) - X 2(t) f dt'Wit'r'xiit'W) , (64) 



with a "source" term 



s(t) = eq(t)xh(t) (65) 



and with the Wronskian 
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W(t) = Xi(*)Xa(«) -X2(t)Xi(t) - |sin(2^) (66) 

(making use of tp = (pi = — ^2)- As mentioned in Section 3, the Wronskian is time-independent. 

Inserting the expressions for the source ( |65|) and for the mode functions xi an d X2 (see (|63|)), and neglecting the 
contribution of the decaying mode \i m the source, we obtain 

Xq {t) = c l£ ^- Vcos(|t + V )h{t) - aeW^e-^cosi^t ~ <p)I 2 (t) , (67) 

where the integrals I\ and I2 are 

I^t) = / dt'q(t')cos(-t' + ip)coa(^t' + <p) (68) 

and 

I 2 (t) - J t dt'e 2 ^q(t')cos 2 (^t' + <p) . (69) 

In estimating the magnitudes of the integrals I\ and I2 we will for the first time make use of our assumption 
of random noise. To be more specific, we will model the noise function n(t) (see (0)) as a random walk with unit 
amplitude and step length T _1 (note that this is a time interval!). By inserting into ( |68| ) and using the standard 
formula for the "radius" of a random walk in terms of the individual random step length, we obtain the estimates 

h - g^{tY) x / 2 Y- 1 (70) 

and 

h ~ S wV' rt (~) 1 / 2 r- 1 (71) 

(where we have set the initial time ti = to simplify the notation) . We conclude that the contribution of I\ dominates 
in ( |67|) as long as we consider time intervals t > fi^ 1 . Thus, we get the following estimate for \ q : 

Xq (t) ~ deg^W-^itT) 1 ^- 1 , (72) 

which must be compared with the amplitude of the "homogeneous" mode 

Xh(t) ~ c ie "*. (73) 

By demanding that the contribution of \q be smaller than Xh over a typical time interval for resonance (i.e. for 
t ~ /i _1 ), and inserting the value ( |66| ) of the Wronskian, we obtain as the condition under which the noise has a 
negligible effect on parametric resonance: 
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or 

g < — M ( M r)- 1 / 2 sin(2^) . (74) 
To discuss the consequences of this condition, recall the expression (|58| ) for fj,. The value of fi is maximal in the center 



of the resonance band and vanishes at the band edges. Hence, we conclude from (74) that any noise will tend to 
slightly decrease the width of the resonance bands. However, resonance at the central value of k is not significantly 
effected unless 

g > i(^) 1/2 shu>), (75) 

which in general yields a value greater than 1. Note, however, that for g > 1 the Born approximation is an invalid 
perturbative expansion. 

We conclude that even large amplitude noise is unlikely to interfere with parametric resonance, and that in fact the 
shorter the time period of the noise, the less likely the noise is to influence the resonant modes (this is reminiscent of 
the Riemann-Lebesgue Lemma). 

VII. RESULTS OF NUMERICAL STUDIES 

At this point we show that the results obtained in the previous sections based on analytical approximate methods 
are perfectly consistent with the numerical analysis of the problem. We also verify that the results of previous sections 
hold for both the narrow resonance and the broad resonance cases, which we studied separately. 

Following the notation of previous sections, we have chosen p(u)t) — Acos(2ujt), where A is a constant, and taken the 
noise q(t) — g A n{t) with g being a positive constant and with nit) being a random function of time with characteristic 
time r _1 and amplitude 1. To simplify the analysis we redefine the time variable to a dimensionless time t — ► t = ujt 
so that the evolution equation for \k reads 

X+ [£ 2 + Acos(2t)+ 5 An(t)]x = 0, (76) 

where A = A/uj 2 and E 2 = ujI/lo 2 . 

In the narrow resonance regime, the first resonance band [E 2 = 1) is the most important. This is because both the 
width of the N-th band and the correspondent Floquet exponent are 0(X N ). This implies that the value of E 2 can 
be shifted (e. g. due to the presence of the noise) by an amount comparable to (but smaller than) A without moving 
a resonant mode into a stability band. Only if the amplitude of the noise function is large (g > 1/2) it is possible 
that the resonance is significantly affected (see Fig. 1). 
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To be explicit let us recall the main characteristics which differentiate the two regimes of parametric resonance 
analysed here: Narrow resonance is defined as the regime where |A| <C 1; E 2 = N 2 + O(X), N = 1, 2, and Floquet 
(Lyapunov) exponent /i <C 1. Broad resonance is characterized by |A| > 1; E 2 > |A and \i ~ 0(1). In the broad 
resonance regime A can be much bigger than 1 but the width of the instability may not be of the same order of A , so 
that even a small change in E 2 (when compared with E 2 itself) can be enough to shift a mode out of the resonance 
band. This implies that even noise with g <C 1 can be important. Moreover, since the growth rate is large (/i = O(l)), 
the particular behavior of %(t) is much more sensitive to the functional form of the noise in the broad resonance than 
in the narrow resonance regime (see Fig. 2 and the discussion below). 

In our numerical work, equation ( |76| ) was solved using a MATLAB integration routine, with a noise function which 
was taken to be piecewise constant in intervals of length and whose amplitude was chosen at random from a 
uniform distribution on the interval (— 1, 1). The most important numerical results can be summarized as follows (all 
analyses were done for g < 1): 
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FIG. 1. Sample of narrow resonance regime results. Xk(t) is plotted (vertical axis) against dimensionless time (horizontal 
axis) for a value of k in the first resonance band, E 2 = 1, and for A = 0.4. 

(a) No noise is present, q(t) = 0. 

(b) The noise is characterized by q ~ 0.010 A, V = 10 and g = 0.8 (large amplitude and large noise frequency). 

(c) The noise is characterized by q « 0.028 A, V = 1/40 and g = 0.4 (small amplitude and small noise frequency). 

(d) The noise is characterized by q w 0.433 A, V = 1/20 and g = 0.9 (large amplitude and small noise frequency). In this case, 
the mode k is shifted out of the resonance band. This is an example for which the noise is important. 
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(i) The noise has practically no effect when T, the inverse of its characteristic time, is much bigger than the frequency 
lu of the periodic driving function p(u)t), i.e. when T 3> u>. This is shown in Figures 1(b) and 2(b) where r » 1, 
and it is consistent with the results of Sections V and VI (see e.g. Eqs. ( f74| ) and ([75])). On the other hand, if T is 
small compared to to, then the ergodic hypothesis is not satisfied, and the noise may completely alter the efficiency of 
energy transfer (see e.g. Figures 1(d) and 2(d)). 

(ii) In the case of narrow resonance, the noise is effective only if its magnitude is of the same order or larger than 
the amplitude of p(u>t), i.e. g > 1/2. For values of g significantly smaller than 1 (g < 0.5, say) the noise is not effective 
for "central" modes in a resonance band. This result was obtained in Sections III, IV and V and is illustrated by 
Figure 1(c). 

(iii) The presence of noise may lead to a much faster exponential growth of a mode initially in a resonance band. 
This effect is particularly important in the case of broad resonance (Figure 2(d)). 

(iv) After having analyzed the time evolution of Xk (t) for many different random walks we verified that an important 
quantity to be considered is q, the end point distance of q(t) divided by the number of steps, over the total time 
interval of reheating t. Note that q — gXh ~ g\/ (tT) 1 ' 2 . If q is small compared to the amplitude A of p(u>t), i. e. if 
g/itYfl 2 < 1, then the resonance is preserved. For some very special random walks in the broad resonance regime 
(see Figure 2(c)) there are specific realizations for which the noise decreases the Lyapunov exponent, but when taking 
the average over many realizations, agreement with the analytical results is restored. 

In order to understand the above results, notice that parametric resonance of a given mode of the field x disappears 
if the noise is able to keep the mode out of the resonance band during a sufficiently large time interval (when compared 
to the total evolution time). 
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FIG. 2. Similar plots as in Figure 1, but now for the broad resonance regime. We consider a mode in the tenth resonance 
band, E 2 = 100, and we choose A = 100. 

(a) No noise is present, q(t) — 0. 

(b) The noise is given by q ps —0.004 A, V = 12.5 and g = 0.2 (i.e. large inverse characteristic time and large noise amplitude). 

(c) The noise is given by q ~ 0.013 A, with the same F and g as in the case (b). 

(d) The noise is given by q ~ —0.025 A, V — 0.2 and g — 0.30 (small F and small noise amplitude). In this example the 
noise has an important effect. It greatly increases the rate of particle production. In this particular case, the noise n takes 
on negative values for a few successive time steps, which keeps the mode in an instability region whose Floquet multiplier is 
larger than without noise. Thus, even though q is small, the overall effect is an increased resonance strength due to the very 
fast exponential growth during a short time interval. 



To see exactly what this means let us write the band number N in terms of the parameters of Eq. (|76|). From 
the theory of the Mathieu equation we have (neglecting the noise) N Q — \E\. According to the point (iv) above, the 
presence of the noise can be approximately taken into account by shifting the value of E 2 to E 2 + gXn = E 2 + q/to 2 . 
Then, by defining iV 2 = N 2 + q/uo 2 where N a is the band number when the noise is neglected (N 2 = E 2 ) we obtain 



AN ee N - N = N ^ = N o9 -^n . 



(77) 



This result means that for large N a , that is to say for modes initially in one of the resonance bands with large band 
number, even a very small mean noise (when compared to the energy of the mode u> 2 ) can shift a given mode through 
several bands and occasionally into a stability band. However, since the same noise is assumed to act independently 
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upon all the modes, it is easy to see that some of the modes initially in a stability band can be shifted to an instability 
band. This effect is important in the broad resonance regime, where the resonance occurs for E 2 ~ A 3> 1. For 
instance, in the case of Fig. 2 we have X/E 2 = 1, so that AiV = 5gn. Then, even for g w 0.1 there are special 
samples of the noise that have a significant effect on the resonance of the boson field (a small effect of the noise means 
AN < 1) . In particular, for Figs. 2(b) and 2(c) we find AN ~ -0.02 and AiV ~ 0.08, respectively. 

On the other hand, for small N a in principle the mean noise could be of the order of lo 2 without jumping to the 
next resonance band. However, this corresponds to the narrow resonance regime A = e <C 1, where the width of the 
instability bands are very narrow and the mode can easily be shifted to a stability band. In fact, for a mode k initially 
in the center of the first band the mean noise q must satisfy \q\/ojl < e/2 in order for not to take the mode out of 
the first resonance band. Notice that e is by definition the ratio A/oj 2 so that the (narrow) resonance is preserved as 
long as \q\ is smaller than one half of the amplitude A of the periodic driving function p(uit) = Acos(2i). In the case 
of Fig. 1 we have X/E 2 — 0.4, so that AiV = OAgn, which implies a small effect of the noise even for g of the order 
of 1. For instance, in the case of Figs. 1(c) and 1(d) we have respectively AiV ~ 0.01 and AiV ~ 0.08. 

According to the general analysis of Sect. II, ergodic noise always increases the rate of particle production. In 
the numerical work, however, we find examples where over a small time interval the rate of particle production is 
decreased as a consequence of the noise. This is manifest in the case of small T (Figures 1(c), 1(d) and 2(d)). This 
apparent contradiction disappears once we realize that for small T and for a finite reheating interval, we are in fact 
not taking the average of q{n 1 1) over many realizations in the sample space f2, but instead choosing only a few special 
noise functions for each calculation. By the ergodic hypothesis used in Section II, we should expect that numerical 
results to agree with the exact results only if the reheating time t is very big compared to the characteristic time 
T = r _1 of the noise. (The time t must be very big, compared to T, to ensure that the noise explores all possible 
realizations in the "noise space".) This condition is grossly violated for the parameters of Figs, lc and Id. There are 
also cases (see Figure 2(c)) when T is large but for certain realization of the noise the Lyapunov exponent decreases. 
In all these cases, however, the mean Lyapunov number over many realizations is larger than the Lyapunov exponent 
in the absence of noise. It would be of interest to further study the dispersion of the Lyapunov exponents for identical 
values of the physical parameters. 



22 



VIII. CONCLUSIONS 



Wc have studied some effects of noise on parametric resonance. Specifically, we considered spatially homogeneous 
noise in the time dependent mass responsible for the parametric resonance instability. Assuming that the noise is 
ergodic we showed that the presence of the noise leads to strict increase in the rate of particle production. We 
demonstrated also that the resonance is rather insensitive to the presence of small noise. Under the assumption that 
the time dependence of the noise can be modeled as a random walk with a characteristic step length, we derived 
estimates for the amplitude of the noise for which it can be shown that the resonance persists. We demonstrated that 
even if the dimensionless amplitude of the noise g is of the order 1, resonance is not affected provided the time step 
r of the noise is sufficiently short. In a subsequent letter we will extend these results to the more interesting 
case of spatially inhomogeneous noise. 

APPENDIX 

This last section gives an outline of the mathematical results that are used in this work. Complete proofs appear 
in the literature citations. 

Proposition 1: When q(t;K) is given by a translation invariant ergodic measure dP{n) on C(R), the limit exists, 

iV— >oo iv 1 

Furthermore, for almost every realization k (with respect to the probability measure dP(n)) the individual limits exist, 

N— too IV 1 

and they are equal to JL. 

The proof of the first statement follows from an argument involving the sequence e(n) — E(log \\IFj =1 $> q (jT, (j — 1)T)\\) 
and its subadditive property, e(n + m) < e(n) +e(m). The proof of the second statement follows from the subadditive 
ergodic theorem, ]l9| . 

Proposition 2: The generalized Floquet exponent (Lyapunov exponent) ^,{q) is continuous with respect to q in the 
topology on C(R) of uniform convergence on compact sets. 

This result and its proof may be found in [jlTj and it follows from Sturm - Liouville theory and the nesting 
property of the Weyl limit circles. 
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The third qualitative result has to do with a monotonicity property of the generalized Floquet exponent in ergodic 
systems. Consider the probability space fl = C(R) with two ergodic invariant measures dPi and dP2- By Proposition 
1 the two associated generalized Floquet exponents [i\ and [ii are constant almost everywhere on the support of 
their respective measures. The following result states that the class of problems (|^) with positive generalized Floquet 
exponent is nondecreasing with respect to the support of the measure. 

Theorem 3: (S. Kotani JT^/j 7/supp {dP\) C supp {dP2) and /ii > 0, then /12 > 0. 

In our case dP\ is supported on the periodic function p{t) and its translates, and dPi will be taken to be the 
description of the statistics of the realizations {p(t) + q(t; k)} C C(R). If q(t) = is a possible realization in the 
support of the probability measure dP(n), then the support of dP2 contains the support of dP\. We are therefore in 
the situation described in Theorem 3. 

The final mathematical result of this article is also the central one to our argument. Consider a probability 
distribution dA on the matrices 6 SL(2,R) (in fact the result applies more generally to SL(n, R)). Let Ga be the 
smallest subgroup of SL(2, R) containing the support of dA. 

Theorem 4: (Furstenberg, JT^/J Suppose that Ga is not compact, and that the action of Ga on the set of lines in 
R 2 has no invariant measure. Then for almost all independent random sequences ^ SL(2, R) with common 

distribution dA, 

lim llog(||n^ 1 * j ||) = A>0 . 

iv— *oo iv 

Furthermore, for given v\^v^ £ R 2 , then 

lim ^log((« l! n^ 1 *^ 2 )) = A 

for almost every sequence 

In our setting, realizations {q(t; k) : < t < T} give rise to transfer matrices \E' g (T, 0) G SL(2,R), and the 
probability measure dP on C(R) induces a measure dA on SL(2,R). We will fulfill the hypotheses of the theorem of 
Furstenberg for dA by demonstrating that under our conditions on supp(dP), then Ga = SL(2, R). This will follow 
if we show that the set of random transfer matrices {4 , g (T, 0) : q £ supp(dP)} contains a small neighborhood of the 
identity in SL(2, R), for this must be contained in the support of the induced measure dA. In order to analyse this set, 
observe that the sequence of successive approximations {^j(t, 0)}°^ gives the Taylor polynomials of the fundamental 
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solution ^>(t, 0;q) = ^ q (t,0) with respect to q G C(0,T) about q(t) = 0. Thus the derivative of the transfer matrix 
^(T, 0; q) with respect to q(t) at zero is given by the first term 

W,0;0) • r(t) = £ dt r{t)%\t) ( ^ q ) *o(t) ■ (78) 

From (||) and the fact that det($ (£)) = 1, this is 

For r(i) G (7(0, T) the expression S q ^(T, 0; 0) ■ r is in the Lie algebra sl(2, R). Taking r(t) G supp(dP) in a small 
neighborhood of zero, if 8 q ^>(T, 0; 0) • r spans sl(2,R), then by the implicit function theorem the solutions ^(T,0;q) 
will indeed fill a neighborhood of I G SL(2, R). From the structure of ([79) the rank of d q ^f(T, 0; 0) • r is at most three, 
spanned by {0f (t), 0i(i)0 2 (t), 0f(i)}, and the question is whether these three components are in all cases linearly 
independent in (7(0, T). 

By direct calculation one verifies that products ip(t) = 4>j(t)4>k(t) of solutions of equation (||) satisfy themselves a 
third order differential equation 

^ = -2{jy +P {t)) + (c 2 +P(t))j t w ■ (si) 

At t = the Wronskian for (ill) is 



0102 02 \ / 1 

A(t) _ = I (01) (0i'0 2 ) (0|) = 10), (82) 

"'" 1 (0?) (0i0 2 ) (01)/ V-2(- 2 +pW) 2 

therefore det(A(t)) = 2 for all t, and the three functions 0f (t), 0i(i)02(i), 0|(i) form a linearly independent set. 
Suppose by way of contradiction that the rank of 6 q ^(T, 0; 0) -r is less than three, so that there is a linear relationship 

C10? + C 2 0l(t)0 2 (*) + C30?,(t) = . 

Taking two derivatives, this implies that (01,02,03)* is a null vector for A(t), contradicting the above assertion of 
independence. 

It is therefore only necessary that span(0 2 (i), 0i(t)0 2 (£), 0|(*)) n supp(AP) contain a small neighborhood of the 
origin in order that 

rank{5 9 *(T, 0; 0) • r : r(i) G supp(dP) n P e (0)} = 3. 

This is surely satisfied for any periodic potential p(t) if B e (0) C supp(ciP), which is our hypothesis. 
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